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We investigate theoretically the long-range electrostatic interactions between a ground-state homonuclear alkali-metal dimer 
and an excited alkali-metal atom taking into account its fine-structure. The interaction involves the combination of first-order 
quadrupole-quadrupole and second-order dipole-dipole effects. Depending on the considered species, the atomic spin-orbit may 
T-H be comparable to the atom-molecule electrostatic energy and to the dimer rotational structure. Here we extend our general 
description in the framework of the second-order degenerate perturbation theory [M. Lepers and O. Dulieu, Eur Phys. J. D, 
2011] to various regimes induced by the magnitude of the atomic spin-orbit. A complex dynamics of the atom-molecule may 
Q^take place at large distances, which may have consequences for the search for an universal model of ultracold inelastic collisions 
Jh as proposed for instance in [Z. Idziaszek and P. S. Julienne, Phys. Rev. Lett. 104, 113202 (2010)]. 



1 Introduction 



c/2 'As illustrated by the present issue of the journal' , the field 
_ O "of cold and ultracold molecules is continuously developing in 
many directions of fundamental and interdisciplinary physics. 
^^|Among many developments, researchers can now create large 
'"O^'samples of ultracold molecules which can undergo elastic, 
i_i .inelastic, and reactive collisions with surrounding ultracold 
■atoms^^. The challenge is at least twofold: the achievement of 
" 'sympathetic cooling of molecules down to quantum degener- 
acy^"^, and the control of elementary chemical reaction at cold 
land ultracold temperatures—^. As pointed out by Julienne (see 
for instance Ref. '■*), low-energy atom-molecule inelastic colli- 
sions can be understood as resulting from the dynamics at large 
\^ distances controlled by long-ranges forces, combined to the dy- 
namics induced by short-range chemical forces. The resulting 
rate can be written as a product involving scattering probabil- 
. . ities in both domains. Due to their complex internal structure, 
^ molecules most often offer many open channels for inelastic or 
reactive processes with other systems so that the related proba- 
bility can be assumed equal to one, so that the entire collision 
is uniquely controlled by well-known long-range forces. Just 
like for ultracold elastic collisions determined by a single pa- 
rameter, namely the scattering length, this approach opens the 
way for elaborating universal models for inelastic and reactive 
atom-molecule collisions, as recently proposed by several au- 
thors^"—. In such models, it is assumed that the inelastic rates 
depend solely on the leading term C„/i?" of the long-range 
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electrostatic interaction between the colliding partners. 

Most of the cases investigated in the papers above involve 
atoms and molecules in their electronic ground state. In a se- 
ries of recent papers we studied the long-range interaction be- 
tween an alkali-metal dimer in its ground state with an excited 
alkali-metal atom, in the perspective of modeling their associ- 
ation into an excited trimer induced by a properly chosen laser 
(photoassociation, or PA) at ultracold energies. In Refs.i^s^S 
(hereafter referred to as Papers I and II, respectively) we 
have characterized the first-order quadrupole-quadrupole and 
second-order dipole-dipole interactions - varying as C^/R^ 
and Cq I (where R is the distance between the colliding part- 
ners) and associated with the operators Vqq and V^^ , respec- 
tively - between a homonuclear alkali-metal dimer in the low- 
est level = of its electronic ground state X^E+, and an 
alkali-metal atom in the first excited state ^P. We have located 
a range of atom-dimer distances with an upper bound _R* in- 
side which both kinds of interaction compete with the rotational 
energy of the dimer Potential energy curves asymptotically 
connected to different dimer rotational levels N are thus cou- 
pled, inducing complex patterns further characterized in Ref.— 
(hereafter referred to as Paper III). The lower bound of this re- 
gion is limited by the so-called LeRoy radius below which 
exchange interaction takes place. 

In this series of three papers we have applied our formalism 
to the interaction between a Cs2 dimer and a Cs atom but it can 
be applied to all alkali-metal combinations. Considering other 
species will modify the limits of the \Rlr\ ^^p] region but not 
the overall aspect of the potential energy curves. This is due 
to the moderate variation of the electrostatic and rotational en- 
ergies from Li to Cs (see Table [T]i. In contrast, the spin-orbit 



1 



2 GENERAL FORMALISM 



splitting of the first excited atomic state n^P varies dramati- 
cally from Li to Cs. Lithium is actually the only species whose 
fine-structure splitting is comparable to the rotational and elec- 
trostatic energies. For all the other species, the spin-orbit split- 
ting is much higher. Following^-', we define two distinct cou- 
pling cases in analogy to Hund's cases: 

• From Na to Cs, the fine-structure splitting is so large that 
the two fine-structure components are not coupled by the 
electrostatic interaction. The spin-orbit Hamiltonian Hso 
is part of the zeroth-order energy ^ = Ex. 1,^=0 + 

EnP + (jiso^ of a given state p, where Ex,va=a is the en- 
ergy of the electronic and vibrational dimer ground state, 
and Enp the energy of the atomic n^P state without fine 
structure. The electrostatic interaction is then obtained by 
diagonalizing at each R the Hamiltonian 



WiiR)=BoN^ + Vgg+V^ 



(2) 



(1) 



with Bq the dimer rotational constant in the electronic and 
vibrational ground state. This situation, denoted as case 
"IC" in Ref. will be illustrated with cesium in the 6^P 
manifold. 

For Li in the 2^P manifold, the two fine-structure com- 
ponents are coupled by the electrostatic interaction, and 
the spin-orbit interaction is included in the perturbation 
Hamiltonian 



T>(2) 



W2{R)^BoN' + Hso + Vyg 

The zeroth-order energy of state p reads Ep . 

Ex,vd=o + E, 



(2) 



np. Due to the very weak variation of the 
electrostatic properties between the 2^Pi/2 ™d 2^P;^^2 
states, we will consider that those properties do not depend 
on the fine-structure level. This situation corresponds to 
case "2A" inRef.^^. 

Table 1 Some atomic properties relevant for the present study: the 
mean square radius {r^p) on the atomic lowest np orbital, its 
spin-orbit splitting AEfs, the rotational constant Bo of the lowest 
vibrational level of the corresponding dimer, and the only nonzero 
tensor component 52 of the dimer quadrupole moment at its 
equilibrium distance (see Paper 1) 





Li 


K 


Cs 


(rip) (a.u.) 


27.1 




62.7 


AEf, (cm-') 


0.335 


57.7 


554 




Lia 


K2 


CS2 


q2 (a.u.) 


10.5 


15.7 


18.6 


Bo (cm-i) 


0.673 


0.0567 


0.0117 



In Section |2] we briefly recall the general formalism that we 
used in our previous work, emphasizing on the modifications 



induced by the presence of the atomic spin-orbit. We derive in 
Section [3] the expressions for the first-order and second-order 
long-range interaction of the atom-molecule system, referring 
to the two cases "IC" and "2 A" above. Then we describe our 
results for the Cs*H-Cs2 case (Section |4]i and for the Li*H-A2 
case, A being an alkali-metal species (Section |5]l. Finally we 
discuss the implications of our results in the perspective of the 
universal model proposed in Ref. '"^ '^ (Section|6ll. 

2 General formalism 

We start from the general form of the electrostatic energy be- 
tween two interacting charge distributions A and B, whose cen- 
ters of mass are separated by the distance R, 



Vel{R) 



1 



+00 L< 

E E j11+La+Lb 



LaLbMWl 



SfA)Ql'UfB). 



(3) 



where L< is the minimum of La and Lb- Eq. (O is the well- 
known expansion on the electric multipoles of A and B. Each 
multipole of order Lx (with X = A^B) is associated with the 
tensor operator i^x), which can be expressed in a coordi- 
nate system whose origin is the center of mass of X 



QlUrx) 



47r 



2Lx 



^ iex 



(4) 



where is the value of each charge i composing X, and 
are the usual spherical harmonics. In Eq. (O, the assumption 
has be made that the quantization axis is the one pointing from 
the center of mass of A to the center of mass of B, hence the 
factor JlaLbM reads 



LaLbM 



(^1)^^ {La + LbV. 
^{La + M)\{La~M)\ 
1 

^{Lb+M)\{Lb- M)\' 



(5) 



In our particular case (Fig. [U, the system A is the dimer, the 
system B is the atom, and the quantization axis Z joins the 
dimer center of mass and the atom. 

In what follows, we want to characterize the first-order 
quadrupole-quadrupole interaction (defined by La — Lb = 2 
in Eq. ^) and the second-order dipole-dipole interaction, (de- 
fined by La ^ Lb ~ 1), in the jj coupling case. Namely, the 
atomic state is characterized by the quantum numbers associ- 
ated with the outermost electron: n the radial quantum number, 
£ the orbital angular momentum, j the total angular momen- 
tum, and Lu its projection on the Z axis, joining the center of 
mass of the dimer and the atom. We shall also use s the spin 
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3 FIRST-ORDER AND SECOND-ORDER LONG-RANGE INTERACTIONS 




Figure 1 (Color online) The two coordinate systems, XaYaZa 
(D-CS) and XYZ (T-CS) defined for the dimer and for the trimer, 
respectively. The Za axis is along the dimer axis, while Z is oriented 
from the center of mass of the dimer towards the atom B. The Y and 
Ya axes coincide and point into the plane of the figure. The 
subsystem A in this figure is the Cs2 molecule, the subsystem B is 
the Cs atom. The T-CS is related to the laboratory coordinate system 
(xyz) by the usual Euler angles (q, j3, 7), not represented here. 



of the electron (s = 1/2), its projection a on the Z axis, as 
well as A, the projection of £. Using these notations, the atomic 
spin-orbit Hamiltonian reads 



Hso = ALs, (6) 
where A — 2AEfs/3, and its matrix elements are diagonal 

Hso) = ^ UU + 1) - i{l + 1) - s{s + 1)) . (7) 

The diatomic molecule is in its ground electronic state 
and vibrational level \vd = 0). The quantum numbers associ- 
ated with the rotation of the dimer nuclei are the angular mo- 
mentum N and its projection m on the Z axis. The total angular 
momentum J, associated with the mutual atom-dimer rotation, 
will be considered in a future work. On the contrary, the pro- 
jection 

rj = TO + w (8) 

of J on the Z axis, will be extensively used, as it is a conserved 
quantity. In comparison to Paper III, the atomic part of the 
electrostatic interaction is the only one which is modified. That 
is why we will focus on it in what follows. 

The atomic multipole moments have now to be expressed in 
the jj coupling case. As the electric-multipole-moment opera- 
tors Qff act on orbital part of the atomic state, we first decom- 



pose the atomicstate, labeled \ juj) on the corresponding decou- 
pled basis {| Act)} 



(9) 



where C^'li^g is a Clebsch-Gordan coefficient. Combining the 



matrix element of Q^^ in the decoupled basis, 
(m^iAiaiigf Inz^Aaaa) 



/ 2^2 + J- /~L \ f^liO ^fiAi 



(10) 



where (ffj) = \^n^i^nn2i2j2) matrix element associ- 

ated with the operator for the outermost electron of the atom, 
and the formula^i 



al3S 



dS 

aaf(j) 



-if+^+'^+f ^(2c+l)(2rf+l) 
a b c 



e f d 



c: 



(11) 



where ^ ^ ^ ^ is a Wigner 6-j symbol, and using the 

fact that ^1 + i has the same parity as £2, we get to the expres- 
sion 

^^+^+^"^+^ v/(2^2 + l)(2j2 + l)(rf2) 



= (-1) 



^l20LQ 



^j2U02LM ■ 



(12) 



In case "2A" defined in the previous section, the radial matrix 
element is independent from ji and 32- 



3 First-order and second-order long- 
range interactions 

The matrix element coiTesponding to the quadrupolar interac- 
tion Vqq is obtained by combining Eq. (fT2l l with L = 2 and 
£1 — £2 for the atomic part, and Eq. (12) of Paper III for the 
dimer part, which yields 



NiuiijiUJi 



9? 



A^2 "12^2^2 



24 (-1)^+-+^/^ y^f^V(24 + l)(2j2 + l) 



£ i 

^ 2 

h 2 



J2 



i?5 



E 

A/=-2 



"-^ Af2 ni 2 2 M "-^ J2 tiJ 2 2 - M 

(2 + M)!(2-M)! ' 



(13) 
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3 FIRST-ORDER AND SECOND-ORDER LONG-RANGE INTERACTIONS 



where is the tensor component of the dimer quadrupole mo- 
ment along its internuclear axis Za- The angular factors of 
Eq. ( fT3b impose strong selection rules: (i) mi +bJi = m2+0J2, 
which means that the quantum number Q, (see Eq. (O) is con- 
served; (ii) Ni = N2 ± 2; and (iii) ji, j2 and L = 2 must 
satisfy the triangle rule. It is important to remark that the latter 
selection rule is not satisfied for ji = j2 = 1/2, which is of 
strong importance for the Cs2+Cs interaction (case "IC")- If 
the cesium atom is in the 6^Pi/2 fine-structure level, it only in- 
teracts with Cs2 through the second-order dipolar interaction, 
which is certainly not favorable for the existence of long-range 
vibrational levels below the 6^Pi/2 dissociation limit of the 
trimer. 

The second-order dipolar interaction is associated with the 

^ (2) 

operator V^^ , whose matrix elements can be written as func- 
tions of the dynamical polarizabilities of the two fragments at 
imaginary frequencies (see Eq. (12) of Paper II) 



N2m2j2i^2 ) 



we get to the final expression 



-M-M 



2 E 



{En 



En 



X (2^ + 1) ^(2ji + 1) (2i2 + 1) 



X {'''nljin'l'j') {rn£j2n'i'j') [^£010 



32 



f 1 e' f\ f 1 f 



(17) 



M 

\ /• + 00 



The key point is now to connect Eq. ( [TtI i to the isotropic po- 
larizability of the atom, which will depend on the considered 
cases "IC" and "2A". 

In the case "IC" illustrated by cesium, the isotropic polar- 
izabilities are different for the two fine-structure levels ^^Pi/2 
and &P'i/2- Moreover, as the subspaces associated to those two 
levels are fully decoupled, the matrix element of V"^^^ are zero 



— / da;aX}M'(*'^)a-M-Af'(*'^) 







if ji 7^ i2- As pointed out in Paper II, the isotropic polarizabil- 
ity oinij corresponding to the level nij with sublevels u) (not to 
be mixed up with the frequency) is obtained from Eq. ([TtT i by 
carrying out a sum over all values of oj' and an average over uj 



AM' 



(14) 



where Q{x) is Heaviside function, and the letter b = n'£' f 
stands for all the quantum states of the atom accessible through 
dipolar transitions. The last two lines of Eq. (fT4l l are contri- 
butions due to the downward atomic transitions. They depend 
on the dimer dynamical polarizability a^l^^ at the (real) fre- 
quencies of the atomic transitions, which is given in Eq. (20) 
of Paper III. 

In Eq. (O, we have introduced the atomic dynamical polar- 
izability in the coupled basis 

= 2(-l)^^V {Eb-EnijJ 



b {Eb - Enlji ) - z 



$6 n£j2Uj2)Al5) 



where z can be either real or complex. The first line of Eq. ( fTSl l 
depends on ji, but not on j2'- the subtle reason for this will be 
explained in the what follows. By applying Eq. ( fT2] ) to L = 1, 
and putting the primes in upper indices of the Clebsch-Gordan 
coefficients using the identity 



anij (z) 



E 



{En'i'j' — Entj ) 
~ n'/j> [En'l']' - Enljf - Z?- 

X (2^+l)(2/ + l)(f„^,W£'j')' 
2 



X C, 



f 1 e' 



(18) 



Here ac is the polarizability of the atomic core which is the 
same as in Paper II, and we used the identity J^a-y (^aab/s) ~ 
2IX1 ■ situation is thus similar to the one of Paper II: be- 
cause of angular factors, the quantity a^'^^jvf cannot be re- 
lated to the sole polarizability of the nij atomic level, and 
contributions from different j j' transitions must be sep- 
arated. For example, the contributions j — 3/2 j' = 1/2, 
j = 3/2 / = 3/2 and j = 3/2 j' ==5/2 for cesium in 
the 6^^3/2 state must be distinguished. Following Eq. (24) of 
Paper II, we express the isotropic polarizability as 



E ^nijn'e'j' + ttc , 
n'i'j' 



(19) 



(-1) 



a+2b-c-fi 



2c+ 1 



2a 



-<^C76-,9 : 



(16) where the state-to-state polarizability Unejn'i'j' is related to 
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4 CASE"1C": CS*+CS2 



-M-M 



_ or 2ji + 1 

' ~ '^"^■^^'^ ^ 27' + 1 

e'={i-i,e+i) n' 

+3' 

LJ2IM' 



X 



(20) 



It is important to keep in mind that the quantum numbers I' 
and j' are related to each other by the condition £' — s < j' < 
£' + s. So, the possible transitions are: — ^ n''^Si/2 and 



and 11? P-- 



3/2 



for an alkali-metal atom in a n P1/2 state, 
Si/2. n^Pi 



3/2 



n''^D3/2 and n2p3/2 



, for an alkali-metal atom in a n^P3/2 state. 
The case "lA" corresponds to a lithium atom in the 2^P 
state. As we have not calculated the dynamical polarizability 
of lithium, we will only present the formalism and no numeri- 
cal results. We assume that the different fine-structure levels of 
the n£ manifold can be coupled by electrostatic interaction, but 
that the energies of the states n£j and n'£'j, , and the transition 
dipole moments from state n£j to state n £'j, do not depend on 
j and j'. It is thus more relevant to express the isotropic polar- 
izability ctne in the Russel-Sanders coupling scheme, which is 
given in Paper II, Eqs. (22) znd (24), 



nl) 



( 



'^eoio 



(21) 



In order to calculate a^j'^i (z), it is, like in Paper II, necessary 
to separate the different £ £' transitions (here P —i' S and 
P D). By introducing the state-to-state polarizability. 



Olnln'l'{z) 



[En' 



E, 



3 {En' 



X [rnin't 



- E„ 



■i'O 

mio 



(22) 



we finally get to the relation 

l' = {l~l,i+l) n' 
J<+1 

X (2^+1) V(2ji + l)(2j2 + l) 



j'=j> -1 



^ 1 n 



+3' 

^jiwilA/^j2tJ2lM' ' (23) 






(b) 




N= 5 . 


/ 1 ^''^ • 





Atom-dimer distance R (a.u.) 

Figure 2 Long-range potential energy curves between a ground-state 
Cs2 and an excited Cs(6^Pj=i/2) atom, as functions of their mutual 
separation i?, for = 1/2 and for: (a) the even values of A'^ and (b) 
the odd values of A'^. The origin of energies is taken at the 
CS2(X, = 0, TV = 0)+Cs(6^Pi/2) dissociation limit. 



with j> the maximum of ji and j2, and j< their minimum. 



4 Case "IC": Cs*+Cs2 

In this case the two atomic fine-structure subspaces are fully 
decoupled and treated separately. Figs.|4]and|5]present the po- 
tential energy curves characterizing the interaction between a 
Cs2 dimer in the lowest rotational levels N and a Cs atom in 
the &P\i2 level, obtained after diagonalization of W\{K) (see 
Eq. ([T]i). The curves are sorted by values of |ri| (here equal to 
1/2 and 3/2) and parity of iV. Table |2] presents the correspond- 
ing Ce coefficients for = and 1, all the C5 coefficients 
being zero. The method used to calculate the polarizabilities of 
Cs2 and Cs are the same as those described in Paper II. 

The potential energy curves shown on Figs. |2] and [3] have 
very similar features. Most dimer rotational levels split into 
two curves as the two fragments get closer to each other. Typ- 
ically, for i? > 80 a.u., one curve, characterized by a negative 
Ce coefficient, is attractive, while the other, characterized by 
a positive Cg coefficient, is repulsive. The positive Ce coeffi- 
cients are due to the highly-negative parallel polarizability of 
the dimer at the atomic &P\i2 6^5 transition frequency 
(-5160 a.u.), compared to -3037 a.u. at the &P-ii2 ^ G^S* tran- 
sition frequency. For lower atom-dimer distances, the repulsive 
curves turn attractive, due to the coupling with the attractive 
curve connected to the higher dissociation limit. The result- 
ing long-range potential barriers, whose height can go up to 
0. 1 cm^^, could prevent collisions in the ultra-cold regime. 
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5 CASE "lA": LF+Aa 




0.4 
0.3 

I 0.1 

LU 



-0.1 











W= 3_ 


■ / / ./ / . . . , w 



100 

Atom-dimer distance R (a.u.) 



Figure 3 Same as Fig. [21 for = 3/2, j = 1/2. 



iV 




(a.u.) 





1/2 


-11022 


1 


1/2 


-19952 


1 


1/2 


6840 


1 


3/2 


-19952 



Table 2 The Ce coefficients of the 

Cs2(X^E+, «d = 0, Af)-l-Cs(6^Pi/2) long-range interaction 
calculated for A'^ = and 1. In analogy to a diatomic molecule, the 
states are sorted by absolute values of the total angular momentum 
projection on the Z axis. All the C5 coefficients are zero. 




/ 




" W=5 




a;=3 




^' "I 




^^"^ . , (")' 



100 

Atom-dimer distance R (a.u.) 



Figure 4 Same as Fig. for = 1/2, j = 3/2. 



For Cs in the 6^^3/2 level, the potential energy curves 
(Figs.|4]and|5j look very different from the previous ones. Like 
those obtained in Paper III without spin-orbit interaction, most 
of them are attractive, and give birth to complex couplings for 
i? < 80 a.u.. Those couplings can consist of the two-curve 
crossings, like the ones described in Paper III, but also of three- 
curve crossings. The most striking example of such a crossing 
is pointed out by an arrow on Fig.|3b), but this feature is quite 
general. They will be described in more details in the next para- 
graph, with Li2+Li as a direct comparison with the curves in 
the underlying spinless symmetries will be possible. 

The corresponding C5 and Cg coefficients are presented in 
Table[3] for = and 1 . The square radius of the 6733/2 orbital 



'3/2 ' 

of the cesium was calculated in our group with a Dirac-Fock 

method, which gives (^1p^^^ ~ 78.45 a.u.. The quadrupole 

moment of Cs2 is the same as in Paper I (q^ = 18. 56a. u.). The 
C5 and Cg coefficients are of the same order of magnitude as 
those obtained in the LS coupling case (see Papers I and II). 
The most attractive curve, found for iV = 1 and = 1/2, 
has inherited its behavior from the most attractive E+ curve 
discussed in Paper I and II, although its C5 coefficient is less 
attractive (-1131 a.u. with respect to -1674 a.u.). The associ- 
ated Cg coefficient is positive, but much lower the one in LS 
coupling case (2500 a.u. compared to 51249 a.u.). By contrast, 
all the other Cg coefficients are negative. 



5 Case "lA": Li +A, 



Now the two atomic fine-structure levels are so close to each 
other that they can be coupled by the electrostatic interaction. 
This is illustrated by a lithium atom in the 2^P state and both 
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5 CASE "lA": LI*+A2 







A/=4- 




/y = 2" 




^ N= 



Atom-dimer distance f? (a.u. 




Figure 5 Same as Fig. [21 for \n\ = 3/2, j = 3/2. On panel (b), the 
arrow points out the three-curve crossings discussed in the text. 



(a) 

A/=4,/=3/2- 




\ \ 


C'). 




_ N=3,y=3/2" 




N=3,y= 1/2 ■ 
^ N=1,y=3/2 ■ 


^=1,;= IS- 



Atom-dimer distance R (a.u.) 

Figure 6 (Color online) Long-range potential energy curves (without 
dipolar interaction) between a ground-state Li2 and an excited 
Li(2^P) atom, as functions of their mutual separation J?, for 

= 1/2 and for: (a) the even values of A'^ and (b) the odd values of 
TV. On panel (a), potential curves calculated by neglecting the fine 
structure of Li(2^P) are also displayed: in the (crosses), E~ (full 
squares) and 11 symmetries (empty squares). The arrow points out the 
three-curve crossings discussed in the text. 



N 



\^\ 



1/2 
3/2 
1/2 
1/2 
1/2 
3/2 
3/2 
5/2 



2 

6P3/2 



gO) Cs (a.u.) Cg (a.u.) 



_6_ 

25 






5 + lj 


(v/5-1) 

9-vT05 



25 
9+VT05 
25 



25 






-50312 





-28580 


-1131 


2500 





-35538 


432 


-72712 


-73 


-39759 


1121 


-67219 


-349 


-23944 



Table 3 The Cs and Cs coefficients of the 
Cs2(X^E^, «d = 0, A'^)+Cs(6^P3/2) long-range interaction 
calculated for = and 1. In analogy to a diatomic molecule, the 
states are sorted by absolute values of the total angular momentum 
projection |f2| on the Z axis. As well as in Paper I, the values of Cs 
are given in scaled units of < 



P3/2 



^2- 



by Li2 and K2 dimers. The results presented in this section 
are obtained by diagonalizing W2 (see Eq. (|2]i), but without the 
second-order dipole-dipole interaction. The square radius of 
the 2j) orbital of the lithium atom was calculated in our group 
with an Hartree-Fock method, which gives = 25.50 a.u.. 
The quadrupole moments of Li2 and K2 are equal to 10.73 and 
15.69 a.u., respectively^^. First, the potential energy curves be- 
tween Li2 and Li are displayed on Fig. |6] for = 1/2. On 
panel (a), for the even values of the rotational quantum number 
N , the curves of symmetries Yi^ and H, calculated without the 
atomic fine structure, are also plotted. As the splitting between 
two subsequent rotational levels of Li2 , at least 6i?o ~ 4 cm^ ^ , 
is much higher that the fine-structure splitting (0.3 cm^^), the 
description in LS coupling case seems adequate, especially in 
the coupling region, for 20 < i? < 40 a.u., where each curve 
belonging to the = 1/2 can be clearly identified with a 
curve belonging either to the or to the H symmetry. We 
can also see three-curve crossings (see arrow on panel (a)), like 
in the case Cs2+Cs. This crossing concerns one E+ and two 
n states. Outside the crossing, one of the E+ states and the 11 
state are almost degenerate, whereas, at the crossing, the two 11 
states avoid each other The mechanism is the same for cesium, 
even if the curves resulting from the calculation in the jj cou- 
pling case do not have such strong components coming from 
one given LS-coupling-case symmetry. 

At last, we consider the interaction between Li(2^P) and 
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6 DISCUSSION 




Atom-dimer distance R (a.u.) 



Figure 7 Long-range potential energy curves (without dipolar 
interaction) between a ground-state K2 and an excited Li(2^P) atom, 
as functions of their mutual separation R, for |n| — 1/2 and the even 
values of A'^. 



K2 for which potential energy curves are plotted on Fig. [T] for 
|ri| = 1/2. This situation is particularly interesting since two 
asymptotic channels namely N = 0,j = 3/2 and N ~ 2,j = 
1/2 respectively at 0.335 cm^^ and 0.340 cm~^, are almost 
degenerate. The two channels can thus be coupled at large dis- 
tances: for instance at R = 120 a.u., the three curves con- 
nected to those two limits and taken with increasing energies 
decompose as (0.2;0.8), (1;0), and (0.8;0.2) on the channels 
(TV = 0, j 3/2; N = 2,j = 1/2). Moreover, the coupling 
region for 30 < i? < 70 a.u. is characterized by numerous 
avoided crossings. 

6 Discussion 

The results displayed in the preceding section clearly show that 
a complicated dynamics is likely to occur at large distances, 
namely well beyond the range of chemical forces, during the 
cold collision between an excited atom and a ground state 
molecule. This situation represent a new possibility to inves- 
tigate the range of validity of the universal model for inelas- 
tic collisions at ultracold energies '^li^. Based on Multichannel 
Quantum Defect Theory (MQDT), this model considers that 
the rate for an ultracold inelastic collision proceeds first from 
the long-range interactions which controls the probability for 
the system to reach the range of chemical forces at short dis- 
tances. A long-range potential with a dominant term C„/i?" is 
characterized by an interaction length a{n) ''^'^^ 



a{n) 



cos 



(a) 



r 



n-l 
n-2 



(24) 



h"^ /(2fj,a{n)'^) is associated to this length. For collision en- 
ergies E = h'^k^/2fi dominated by s-wave scattering, i.e. 
ka{n) << 1, the inelastic rate is independent of k and writes 
K™{n) — 2ha{n)/ii^*. In the present case, the dominant 
term for the interaction between a ground state Cs2 molecule 
and a Cs(P3/2) atom corresponds to n = 5, and for a typical 
value C5 = —1131 a.u. (see Table [3]), one obtains a(5) = 
0.364506(2AiC5/;i^)i/3 ^ 260 a.u.. For a Cs(Pi/2) atom 

(n = 6) one has a(6) = 0.477989(2^C6/?i^)^/* ~ 103 a.u.. 
In both case, the characteristic length is larger than the inter- 
mediate range of distances identified above where spin-orbit 
coupling, rotational energy and electrostatic energy all compete 
together Therefore such an ultracold inelastic collision will be 
actually defined by three domains: (i) the large distances where 
the sole C„/i?" term controls the scattering, (ii) the intermedi- 
ate distances above with a specific treatment of the dynamics, 
(iii) the short distances where one can reasonably assume that 
so many channels are open that the reaction probability is equal 
to unity. The dynamics in the intermediate range could well be 
treated within a coupled-channel framework in one dimension 
as implemented in Ref. * for ultracold ground state Rb (or Cs) 
atoms and RbCs molecules. Thus it is most likely that the va- 
lidity of the universal model of Refs. '"^-'^ will be limited in the 
present case, just like it has been argued for ultracold collisions 
between ground state KRb molecules in their lowest rovibra- 
tional levels 

Another remarkable expected feature of such collisions will 
be the existence of numerous Feshbach resonances induced 
by bound levels of the trimer close to one dissociation limit 
Cs*H-Cs2(iV) interacting with the continuum related to a disso- 
ciation limit with a smaller value of N. Such resonances may 
enhance the photoassociation probability and the decay down 
to stable molecules through the i?,-transfer of the probability 
density of the system to smaller distances, as it is well known 
for atom-atom photoassociation 

Finally, it is also worthwhile to mention that the energy spac- 
ing between molecular rotational levels is -at least for Cs2- of 
the same order of magnitude than the hyperfine splitting of the 
Cs ground state, which could again induce a complex resonant 
dynamics at large distances, which will the subject of a further 
study. 
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